Phenomenology of High Energy Neutrinos in Low-Scale Quantum Gravity Models 
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We show that neutrino telescopes, optimized for detecting neutrinos of TeV to PeV energy, can reveal 
threshold effects associated with TeV-scale gravity. The signature is an increase with energy of the 
cross section beyond what is predicted by the Standard Model. The advantage of the method is that 
the neutrino cross section is measured in an energy region where i) the models are characteristically 
distinguishable and ii) the Standard Model neutrino cross section can be reliably calculated so that 
any deviation can be conclusively identified. 
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Motivated by the absence of a self-consistent theory of 
quantum gravity and the unresolved hierarchy problem 
between the electroweak scale (10 2 GeV) and the Planck 
scale (10 19 GeV), a great deal of attention has been given 
to theories of low-scale quantum gravity which envision 
significant quantum gravity effects at an energy scale of 
the order of M s ~ 1 TeV WW. In these scenarios, po- 
tentially large effects on high energy processes may occur 
due to the contributions from, e.g., Kaluza-Klein excita- 
tions of gravitons (KK) or other stringy states near M s . 
An interesting motivation for models in which cross sec- 
tions at TeV scale become enhanced is the ultra-high en- 
ergy cosmic ray problem. Protons above the GZK cutoff 
(~ 10 19 eV) interact with the cosmic microwave back- 
ground cataclysmically by the A-resonance ||0]. Thus, 
the cosmic ray events observed above this energy must be 
produced by local sources, or involve new physics. Local 
sources of particles of such energy being unlikely, many 
exotic solutions have been proposed || . A solution which 
has received a great deal of attention in recent literature 
proposes that neutrinos with enhanced cross sections at 
GZK energies constitute the highest energy cosmic rays 
P-p|. This solution requires neutrino- nucleon cross sec- 
tions on the scale of 10's of mbarnes. Unfortunately, most 
scenarios of low-scale quantum gravity as low-energy ef- 
fective theories are valid only up to the order of < M s . 
Above this scale, the naive calculations typically violate 
unitarity g,0|. One has to introduce some ad hoc unita- 
rization scheme, since the fundamental theory, such as a 
realistic string theory, is yet unavailable. It is also very 
difficult to reliably predict the parton distribution func- 
tions needed at GZK energies in neutrino-nucleon inter- 
actions. For these reasons, studies of ultra-high energy 
(~ 10 20 eV) quantum gravity enhancements to neutrino- 
nucleon interactions are extremely speculative. 

These problems are far more manageable at energies 
below or near M s . Unitarity may not be violated at this 
scale, calculations are generally perturbative and the rel- 
evant parton distributions are known at these energies H . 
The characteristics for different theoretical models can be 
also qualitatively distinguishable near and slightly above 



the threshold. Therefore, the TeV regime provides a nat- 
ural scale for probing the features of low-scale quantum 
gravity models. These tests include direct searches in 
colliders such as the Fermilab Tevatron and the CERN 
Large Hadron Collider (LHC) |§-|l|. This letter dis- 
cusses another class of experiments capable of testing 
features of low-scale quantum gravity: multi-TeV to PeV 
neutrino astrophysics. 

Wc have recently witnessed first light of neutrino tele- 
scopes optimized to detect neutrinos in the TeV to PeV 
energy range |13],|l4| . This is the range of laboratory en- 
ergies where the onset of TeV-scale gravity effects on the 
neutrino cross section will first manifest itself. For a neu- 
trino flux $, the number of neutral currents events N v 
observed as hadronic showers in a neutrino detector of 
effective area A is given by the convolution over energy 
of the quantity A 3 / 2 x $ x n x o v . Here n is the density 
of the target that interacts with a neutrino with cross 
section a u to produce a hadronic shower. In our discus- 
sion the detected neutrino flux plays a secondary role. It 
may represent the atmospheric neutrino flux or the flux 
of neutrinos of hundreds of TeV anticipated from gamma 
ray bursts. 

An increase with energy of the cross section for neutri- 
nos to interact with matter beyond a level calculated in 
the Standard Model will signal the onset of new physics 
including the increase anticipated as a consequence of 
TeV-scale gravity effects. It is important to recognize 
that the Standard Model cross section is computed from 
nucleon structure functions probed by HERA experi- 
ments in this energy range. The Standard Model baseline 
against which to measure new physics is known. The 
important observation here is that due to the geome- 
try of the Earth and detector, the angular distributions 
of events with new physics can be significantly differ- 
ent from the Standard Model prediction. These mea- 
surements, though challenging for existing instruments 
like AMANDA and the detectors in the Mediterranean, 
should be feasible for second-generation detectors such 
as IceCube. Unlike first generation neutrino telescopes, 
IceCube can separate interesting high energy events from 



the large background of lower energy atmospheric neutri- 
nos by energy measurement. The instrument can identify 
high energy neutrinos over Air solid angle, and not just in 
the lower hemisphere where they are identified by their 
penetration of the Earth, as is the case with AMANDA. 
Furthermore, due to the enhanced neutrino cross section 
from new physics, the Earth eventually becomes opaque 
to those high energy neutrinos. This leads to the strik- 
ing phenomena that only down-going events from new 
physics at high energies are enhanced. 

For the sake of illustration, we will examine three 
classes of low-scale quantum gravity/stringy models: 

1). ADD scenario [[[J: Large extra dimensions (R < 
0.1 mm) with a flat Minkowski metric. The large 
effects arc due to the high degeneracy of the light 
KK gravitons of mass m KK ~ l/R. The theory 
saturates perturbative unitarity at y/s ^ M a and 
some unitarization scheme has to be introduced. 

2). RS scenario M: One anti-de Sitter dimension with 
a non-factorizable "warped" geometry. The phys- 
ical consequences relevant to our interests are the 
contributions from the KK gravitons of a TeV mass 

& 

3). Veneziano amplitude: Due to the lack of a funda- 
mental theory of quantum gravity, a reasonable pa- 
rameterization to the new physics at the scale M s 
perhaps is to include a sum of possible "stringy" 
states Q. The Veneziano amplitude serves for 
this purpose in describing the stringy physics, be- 
fore non-perturbative quantum gravity effects be- 
come overall dominant. It also manifestly preserves 
unitarity m. 

There are several free parameters in these models. The 
choices for these parameters were selected to illustrate a 
variety of phenomenological features and are not inclu- 
sive. First of all, we take the quantum gravity or string 
scale M s as 1 TeV throughout our studies. In addition, 
the ADD scenario is subject to the number of extra di- 
mensions, although the effect under consideration is not 
sensitive to it |Hj] . The RS model varies with the scale of 
the theory (A = e~ kR M p i) and the mass of graviton res- 
onances (m g ). In our analysis, we will only keep the first 
graviton resonance and ignore the heavier states. Models 
calculated with Veneziano amplitudes are characterized 
by two constants a and b, which parameterize the Chan- 
Paton traces for string models g , in addition to M s . 

Perturbative calculations for the ADD scenario lead to 
a parton-level cross section cr(vg — > vg) ~ s 3 /Mf, which 
violates partial wave unitarity |i?e(a^)| < 1/2 at energies 
above 1.5M S . To have a sensible estimate, partial wave 
amplitudes are taken to saturate unitarity when the uni- 
tarity bound is reached. Calculations are then made with 
only S-wave, with 5 leading partial waves, as well as with 



all partial waves that need to be unitarized. This variety 
of choices reflects our ignorance of how nature chooses to 
restore unitarity at a higher energy scale. However, we 
consider the S-wave unitarization a conservative estima- 
tion for the ADD scenario. The RS model is less likely 
to violate unitary as the rapid growth of amplitudes only 
occurs way above graviton mass scale. The Veneziano 
amplitudes automatically respect unitarity bounds. 
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FIG. 1. Neutrino- nucleon cross sections in a variety of 
models compared to the Standard Model neutral current pre- 
diction. ADD (large extra dimension) models are for all, 5 
and 1 partial waves, up to unitarity saturation (top to bot- 
tom). RS models are shown for A=l TeV, m 9 =500 GeV; A=l 
TeV, m g = l TeV; A=3 TeV, m 9 =500 GeV (top to bottom). 
Models using Veneziano amplitudes are for a = b = 5 and 
a — b = (top to bottom). M s —1 TeV for all models. 

Figure |l| depicts the sum of the neutral current neu- 
trino nucleon cross section within the Standard Model 
and the corresponding cross section due to new physics 
for the different models discussed above. For the rea- 
sons explained earlier, we wish to explore new physics 
not too far from the threshold, and thus do not address 
the energy range above E v — 5 x 10 4 TeV. The in- 
crease in the cross section beyond the Standard Model 
due to TeV-scale quantum gravity starts at a thresh- 
old energy E" ~ 10 3 TeV, corresponding to a neutrino- 
nucleon center-of-mass energy near 1 TeV. Indeed, the 
three models under consideration have distinctive char- 
acteristics near the threshold, as clearly seen in Fig. nl 
For the large extra dimension scenario (ADD), the three 
(solid) curves correspond to those of unitarized ampli- 
tudes (S-wave only, 5 partial waves, and all partial waves 
that saturate unitarity). Up to an energy E u = 5 x 10 3 
TeV, the three solid curves are identical, indicating the S'- 
wave dominance. At higher energies, there will be more 
partial waves to reach the unitarity bound. We consider 
the S'-wave unitarization scheme a conservative represen- 
tation for this model. We see that above E% the effects 
of the new physics in ADD should be clearly visible for 
which the increase in the cross section beyond the Stan- 



dard Model expectations could be orders of magnitude 
at E v = 5 x 10 4 TeV. For the RS scenario, the influence 
with a single KK graviton exchange should be also visible 
if the graviton is not too heavy. The result with A = 1 
TeV and a graviton mass m g — 500 GeV could lead to 
a maximum increase in the cross section of a factor of 
30 at E v = 5 X 10 4 TeV. On the other hand, the effects 
would be unobservable if we take m g — 1 TeV and A = 3 
TeV. The stringy model of Veneziano amplitude M pre- 
dicts an increase about a factor of 6 with a = b = 5 at 
E v = 5 X 10 4 TeV, while the enhancement is negligibly 
small if we take a = b = 0. 
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FIG. 2. Energy distribution of v^ + v^ + v e + v e in- 
duced shower events in IceCube. The upper panels show the 
up-going neutrino events and the lower panels the down-going 
events. The panels are labeled with the corresponding the- 
oretically predicted neutrino flux (see text). In each panel, 



the solid line is the event rate for a 01 '"' f a (S-wave); the 



dashed 3-dotted for <f~ 



+ a 



SM 
Veneziano 



dashed line for a SM + a KS with m g = 
TeV and, the dotted for a SM alone. M, 



with a — b — 5; the 
500 GeV and A = 1 
1 TeV for all models. 



Figures 2 and 3 illustrate the qualitative behavior 
explained above in the IceCube detector. We explore 
three different theoretical predictions of the neutrino 
flux. These have been chosen mostly for illustrative pur- 
poses. The panel labeled charm refers to maximal pre- 
dictions of neutrinos from the decay of charmed parti- 
cles produced by cosmic ray interactions in the atmo- 
sphere [[L5[. The W-B refers to the Waxman and Bah- 
call (W-B) limit on the neutrino flux from astrophysi- 
cal sources that are optically thin to proton-photon and 
proton-proton interactions. This represents a flux of 
£ji = 2x 10~ 8 GeV (cm 2 s sr) -1 g6J. GRB labels 
the neutrino flux from Gamma Ray Bursts, energetic ex- 
plosions in the Universe which occur at a rate of ~ 1000 
per year. The flux accounts for fluctuations in the dis- 
tance to individual GRBs and in the energy they release 
in the form of gamma rays fljllg] . Note that the unusual 
features in the angular distribution of GRB neutrinos are 
a result of these fluctuations |17J] . Although these results 
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FIG. 3. cos(# zon ith) distribution of u^ + P^ + Ue + Ue induced 
shower events in IceCube. cos(0zenith) = —1 corresponds to 
vertical up-going neutrinos, and cos(# ZC nith) = to neutrinos 
coming from the horizon. The panels are labeled with the 
corresponding theoretically predicted neutrino flux that has 
been used to obtain the event rate (see text). The legend for 
different models is the same as in Fig. 2. 

The upper panels of Fig. 2 show the up-going neu- 
trino induced shower events (—1 < cos6* zon i t h < 0) 
and the energy spectrum for down-going shower events 
(0 < cos6 zen ith < 1) is shown in the lower panels. The 
event rate peaks at the neutrino energies at which the 
product $>(E V ) x <j u {E v ) maximizes. Figure 3 shows the 
zenith angle distribution of the neutrino induced shower 
events for the same three neutrino fluxes. Notice that 
the shower rate has a contribution from Standard Model 
charged current neutrino interactions which is present 
in all plots. The angular distribution of the up-going 
events is dominated by absorption in the Earth whereas 
for down-going events absorption is not important and 
the angular distribution reflects the detector's effective 
volume to contained showers. 

Little sensitivity to new physics is anticipated when 
looking for up-going events, while increase in the event 
rate for down-going due to new physics may be visible 
if the neutrino fluxes are large enough at high energies. 
Within the large extra dimension model (ADD), the sig- 
nal is observable even for the conservative scenario with 
the S-wave unitarization. Reasonable sensitivity is ex- 
pected to the RS model with a single graviton exchange 
if TO g < 500 GeV, and similar sensitivity is predicted 
for the Veneziano scenario if a, 6 ~ 5. The excess of 
events due to new physics in the charm case is very 
small (< 1 per km 2 and year) because the flux drops 
quickly in the energy range where the cross section in- 
creases, although there is large uncertainty for the atmo- 
spheric neutrino flux at this energy regime. The W-B 
flux gives a sizeable signal due to the larger neutrino flux 
rate near the threshold energies of new physics. The inte- 



grated numbers of events are summarized in Table I for 
the different neutrino fluxes and new physics scenarios 
studied. The numbers in parenthesis are the event rates 
with muons. They include contributions from muons in 
charged current v^ + £y, interactions within the Standard 
Model and muons produced in hadronic showers induced 
by Ufj, + v^ + v e + v e . The IceCube detector is sensitive to 
down-going neutrino induced muons above 1 PeV since it 
can determine the energy of the events and hence sepa- 
rate them from the background produced by cosmic rays 
in the atmosphere, which is smaller than the neutrino 
induced muon rate for energies above ~ 1— 10 PeV p9[. 
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85 (399) 


10.3 (315) 
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5 (20) 


Large extra-D (all waves) 
Large extra-D (5 waves) 
Large extra-D (1 wave) 


85 (399) 
85 (399) 
85 (399) 


219 (433) 
120 (325) 
112 (319) 


5 (8) 
5 (8) 
5 (8) 


54 (71) 
12 (24) 
9 (22) 


R-S, m g = 500 GeV, A = 1 TeV 

R-S, m g = 1 TeV, A = 1 TeV 
R-S, m e = 500 GeV, A = 3 TeV 


85 (399) 
85 (399) 
85 (399) 


106 (318) 
103 (316) 
103 (315) 


5 (8) 
5 (8) 
5 (8) 


7(21) 
6 (20) 
5 (20) 
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85 (399) 
85 (399) 


103 (315) 
106 (316) 


5 (8) 
5 (8) 


5 (20) 
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FIG. 4. Vf± + z> M + u e + u e shower event rate and fj, event 
rate (in parenthesis) per km 2 yr in IceCube for different the- 
oretically predicted v fluxes and ^-nucleon cross sections, fj, 
events include contributions from muons in charged current 
Vfi + v^ interactions within the Standard Model and muons 
produced in hadronic showers induced by f M + Vy, + v B + v & . 
The muon energy threshold is £,, = 100 GeV, the shower en- 
ergy threshold is 1 TeV and the maximum neutrino energy is 
5 x 10 4 TeV. 

In summary, TeV scale quantum gravity or stringy 
models have distinctive characteristics near the new 
physics threshold and may be probed in large high en- 
ergy neutrino telescopes. The ADD scenario with TeV 
scale quantum gravity has distinctive phenomenological 
features for high energy neutrino scattering. We illus- 
trated this point by a conservative unitarized S'-wave 
approximation. The RS model with a single graviton 
exchange can also provide observable signatures if the 
graviton mass is lighter than about 1 TeV. TeV Stringy 
models using Veneziano amplitudes can produce inter- 
esting features in the case with large Chan-Paton trace 
factors. Studies of these features may allow for discovery 
of such models, or for stronger constraints on the scale of 
new physics. Neutrino astrophysics can provide a means 
to complement searches for TeV scale quantum gravity 
in collider experiments. 
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